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Abstract. The Lyon-Meudon Extragalactic database 
(LEDA) gives a free access to the main astrophysical pa- 
rameters for more than 100,000 galaxies. The most com- 
mon names are compiled allowing users to recover quickly 
any galaxy. All these measured astrophysical parame- 
ters are first reduced to a common system according to 
well defined reduction formulae leading to mean homo- 
geneized parameters. Further, these parameters are also 
transformed into corrected parameters from widely ac- 
cepted models. For instance, raw 21-cm line widths are 
transformed into mean standard widths after correction 
for instrumental effect and then into maximum velocity ro- 
tation properly corrected for inclination and non-circular 
velocity. This paper presents the reduction formulae for 
each parameter: coordinates, morphological type and lu- 
minosity class, diameter and axis ratio, apparent magni- 
tude (UBV, IR, HI) and colors, maximum velocity rota- 
tion and central velocity dispersion, radial velocity, mean 
surface brightness, distance modulus and absolute magni- 
tude, and group membership. For each of these parame- 
ters intermediate quantities are given: galactic extinction, 
inclination, K-correction etc.. 

All these parameters are available from direct connex- 
ion to LEDA and distributed on a standard CD-ROM 
(PGC-ROM 1996) by the Observatoire de Lyon via the 
CNRS. 
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1. Introduction 

This paper gives a detailed description of the reduc- 
tion of astrophysical parameters available through LEDA 
database for more than 100,000 galaxies. It is often re- 
quired by users of LEDA who need a reference where the 
description of parameters reduction is given. Most of these 
reduction procedures were described in previous studies: 

— Central velocity dispersion (Davoust et al., 1985) 

— Kinematical distance modulus (Bottinelli et al., 1986) 

— HI data (HI velocity, flux and 21-cm line width) (Bot- 
tinelli et al., 1990) 

— Diameters (Paturel et al. 1991) 

— Names of galaxies (Paturel et al., 1991a) 

— Group membership (Garcia, 1993) 

— Apparent magnitudes (Paturel et al. 1994) 

— Correction for inclination effect (Bottinelli et al., 1995) 

Here, we will summarize these reduction and give re- 
ductions for some additional parameters: morphological 
types and mean effective surface brightness. 

Each forthcoming section will be devoted to a given 
class of parameter. 

2. General features 

In the following sections, each parameter will be desig- 
nated in a unique way by the word used in the LEDA 
query language. The correspondence between this designa- 
tion and the meaning is given in Appendix B together with 
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the variable name used in programs and the format (FOR- 
TRAN convention). Such a designation aims at avoiding 
subscript or exponant characters which can not be used 
in simple text printing or keyboard entering. 

The error on each parameter is now calculated using 
a method which gives better description of the actual ac- 
curacy of the measurement. In previous catalogs the stan- 
dard error on an average parameter was simply calculated 
from the sum of weights of each individual parameter (the 
weight being the inverse square of individual standard er- 
ror). Now, the standard error is augmented quadratically, 
by the external standard deviation between each parame- 
ter. This estimate will be designated as the actual uncer- 
tainty. The precise expression is explained in Appendix 
A. The main advantage of this new definition is that a 
parameter with low actual uncertainty cannot result from 
discrepant individual measurements. So, this definition al- 
lows the user to select undoubtly good data. n 

All parameters are available through LEDA a or from 
a CD-ROM distributed by the Observatoire de Lyon via 
the CNRS. 

3. Name of galaxies 

We collected the most common names among 40. The dif- 
ferent acronyms used are listed in Table 1 with their ab- 
breviation and the number of occurrences. Some names 
designate several objects (e.g. UGC1 designates two galax- 
ies, PGC00177 and PGC00178). In such cases the name is 
given to both objects between parentheses. 

Since our first catalog of Principal Galaxies (PGC Pa- 
turel et al., 1989 , 1989a) we have added many new galax- 
ies in LEDA database. Each galaxy created in LEDA 
database receives a permanent LEDA number, (with the 
acronym LEDA). Note that LEDA number is identical to 
PGC number for running number less than 73198. PGC 
numbers are sorted according to right ascension and dec- 
lination for epoch 2000. 

Many galaxies are known by their lexical name. 
These names are useful for some nearby large galaxies 
(Dwingeloo 1 and 2; Maffei 1 and 2 etc.). The equiva- 
lence of these names is given in Table 2. 

4. Coordinates 

Equatorial coordinates (Right ascension and Declination) 
are given for two equinoxes 1950 (Besselian coordinates 
aZ1950, del950) and 2000 (Julian coordinates aZ2000, 
cfe2000). Most of published coordinates are B1950 coor- 
dinates. Conversion to J2000 has been made according 
to the "Merits Standards" published in the U.S. Naval 
Observatory Circular. (1983). For computer use, coordi- 
nates are expressed as decimal values (hours to 0.00001 



Acronym 


IN 


Keterence 


PGC 


101258 


Paturel et al., 1989 


MCG 


30662 


T. 7 i T T 1 ' j 1 -1 /"\ /TV -\ r\<—> A 

Vorontsov-Velyammov et al., 1962-1974, 


CGCG 


29825 


Zwicky et al., 1961-1968 


ESO 


17277 


T i j_ 1 Aon 

Lauberts, 1982 


UGC 


13084 


Nilson, 1973 


IRAb 


lli>6o 


IKAb, Fomt bource Catalogue, lyoo 


KUG 


7942 


T~f~\ 1 TV T* 1 * T T 1 /~\ O A 1 /"\/~VO 

Takase & Miyaucm-Isobe ,1984-1993 


SAIT 


7044 


baito et al., 1990 


NGC 


6517 


Dreyer, 1888 


DRCG 


5725 


Dressier, 1980 


IC 


3509 


Dreyer 1895,1910 


T GC 


2573 


jy i j „j_ i i nno 

Karacnentsev et al., 1993 


VCC 


2097 


Bmggeh et al., 1985 


FGCE 


1881 


T 7~ 1 J J 1 1 AAO 

Karacnentsev et al., 1993 


TV /T A n 1S~ 

MARK 


1514 


TV X 1 j l -i r\r>i-f -1 A0 1 

Markaryan et al., 1967-1981 


KCPG 


1206 


Karacnentsev, 1987 


A1N (JIN 


1179 


de Vaucouleurs et al., 197b 


FAIR 


1185 


T71 11 1 C\'—J'—J 1 AOO 

lairall, 1977-1988 


vv 


1164 


Vorontsov-Velyammov, 1977 


nZW 


2714 


Zwicky, 1971 


KARA 


1051 


Karachentseva, 1973 


UM 


652 


Kojoian et al., 1982 


T. HTT r~7T T r 

VIIIZW 


645 


Zwicky et al. , 1975 


ARAK 


595 


Kojoian et al., 1981 


KAZA 


581 


"XT' * -1 A^TA -i r\ 

Kazarian, 1979-1983 


DCL 


570 


Dickens et al., 1986 


Ann 
ARr 


561 


Arp, 1966 


HICK 


464 


TT' 1 -1 AAO 

Hickson, 1993 


UGCA 


441 


Nilson, 1974 


FCC 


340 


T71 001 1 r\f\r\ 

lerguson & bandage, 1990 




291 


Karachentsev et al., 1993 


SBS 


284 


Markanan, 1983-1984 


DDO 


242 


1 lsher & lully, 1975 


tt, TTT'TIVT 

WE1JN 


207 


T.T, 7 ■ 1. 1 aoa 

Weinberger, 1980 


TOLO 


111 


Smith et al., 1976 


SXD 


/ 


rvoou 06 ijauin, lyui 


nSZW 


58 


Rodgers et al., 1978 


MESS 


40 


Messier, 1781 


POX 


24 


Kunth et al., 1981 



telnet lrnc.univ-lyonl.fr - login: leda 



or: http://www-obs.univ- lyonl.fr/base/home_base.fr. html 



Table 1. List of acronyms 



and degrees to 0.0001 for Right Ascension and Declina- 
tion, respectively) . The standard deviation of coordinates 
is generally not known. Thus we are using a flag ipad to 
tell if the standard deviation is smaller than 10 arcsec or 
not. We collected systematically accurate coordinates in 
literature (see Paturel et al. 1989). Recently, we added ac- 
curate coordinates directly obtained from images stored in 
LEDA (Paturel et al., 1996) and from COSMOS database 
(Rousseau et al., 1996). Among the 100872 galaxies 69165 
have accurate coordinates (69 percent). 

Galactic coordinates 12, b2 are calculated (in de- 
grees to 0.01 deg) from all950 and del950 using the 
coordinates of the galactic pole aZ1950(pole)=12. 81667 
(iel950(pole)=27.4000 and the coordinates of the ori- 
gin aZ1950(origin)=17.70667cfel950 (origin)=-28.9167 ac- 
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Lexical name 


Usual name 


PGC number 


LMC 


ESO 56-115 


PGC 


0017223 


SMC 


NGC 292 


PGC 


0003085 


Maffeil 


UGCA 34 


PGC 


0009892 


Maffei2 


UGCA 39 


PGC 


0010217 


Circinus 


ESO 97- 13 


PGC 


0050779 


Sextans A 


MCG -1-26- 30 


PGC 


0029653 


SextansB 


UGC 5373 


PGC 


0028913 


Carina 


ESO 206- 20A 


PGC 


0019441 


Draco 


UGC 10822 


PGC 


0060095 


Fornax 


ESO 356- 4 


PGC 


0010093 


Sculptor 


ESO 351- 30 


PGC 


0003589 


UrsaMinor 


UGC 9749 


PGC 


0054074 


Phoenix 


ESO 245- 7 


PGC 


0006830 


LeoA 


UGC 5364 


PGC 


0028868 


Pegasus 


UGC 12613 


PGC 


0071538 


WLM 


MCG -3- 1- 15 


PGC 


0000143 


Malinl 




PGC 


0042102 


HydraA 


MCG -2-24- 7 


PGC 


0026269 


CygnusA 


A/Tf^r 1 7/11 Q 
Ivl^'Li f -41- O 


PGC 


0063932 


HerculesA 


MCG 1-43- 6 


PGC 


0059117 


Dwingelool 




LEDA0100170 


Dwingeloo2 




LEDA0101304 



Table 2. Galaxies known by their lexical name 



cording to Blaauw et al. (1960). These galactic coordi- 
nates are used to estimate the galactic extinction ag con- 
verted to Burstein-Heiles system (Burstein and Heiles, 
1984) from the relationship given in the Second Refer- 
ence Catalog (de Vaucouleurs, de Vaucouleurs and Cor- 
win, 1976; p32, rel. 22; hereafter RC2). In fact, for galactic 
latitude 62 > 20 deg the difference between both systems 
is negligibly small (except for the zero point difference 
of 0.20 mag due to the fact that Burstein-Heiles give no 
absorption at the galactic pole). In Fig. |l| the difference 
between ag from RC2 and from Burstein-Heiles is plotted 
vs. the galactic latitude. 

The conversion from RC2- to BH-system is the follow- 
ing: 

if | 62 1 < 20 deg: 

ag{BH) = ag{RC2) + 0.70 - 0.045 \b2\ (1) 
and if |62| > 20 deg: 

ag(BH) = ag(RC2) + 0.001 \b2\ - 0.22 (2) 

The galactic extinction is higher than the one pre- 
dicted by RC2 formula for low galactic latitude (62 < 
20 deg). The adopted galactic extinction is ag = ag(BH) + 
0.20, where ag(BH) is calculated from Rel. 1 and Rel. 2. 

Supergalactic coordinates sgl, sgb are calculated (in 
degrees to 0.01 deg) from 12, 62 using the coordinates of 
the supergalactic pole /2(pole)=47.37 deg 62(pole)=6.32 
deg and the coordinates of the origin Z2(origin)=137.37 
deg; 62(origin)=0 deg according to de Vaucouleurs et al. 
(1976). 




20 40 60 8C 

absolute value of galactic latitude |b2 



Fig. 1. Difference of galactic extinction from RC2 and 
Burstein-Heiles vs. the absolute value of the galactic latitude. 

5. Morphological type, luminosity class, and lumi- 
nosity index 

Morphological types (E, SO, Sa ... Sm, Irr) have been 
entered in LEDA as an internal numerical code. This code 
will be treated as a continuous quantity. 

Obviously, the definition of what is a given morpho- 
logical type is not the same for each astronomer. We thus 
adopted the RC3 type code system as a reference one and 
converted to it all type codes of others catalogs. A rms 
dispersion can be attached to each type code for each ref- 
erence allowing us to calculate a weighted mean morpho- 
logical type code t and its actual uncertainty st. In addi- 
tion, some features have been coded in LEDA: ring, bar, 
interaction (or multiplicity), and compactness. These fea- 
tures are simply given as a flag R,B,M for the first three 
parameters and C or D for the last one (C for compact 
and D for diffuse). Further, the numerical code t and the 
features above are used to produce a literal Hubble type 
typ (e.g. SBa). The ranges of definition are in Table 4. 



range of t 


typ 


range of typ 


type 


-5 < t < -3.5 


E 


3.5 < t < 4.5 


Sbc 


-3.5 < t < -2.5 


E-SO 


4.5 < t < 6.5 


Sc 


-2.5 < t < -1.5 


SO 


6.5 < t < 7.5 


Scd 


-1.5 < t < 0.5 


SOa 


7.5 < t < 8.5 


Sd 


0.5 < t < 1.5 


Sa 


8.5 < t < 9.5 


Sm 


1.5 < t < 2.5 


Sab 


9.5 < t < 10 


Irr 


2.5 < t < 3.5 


Sb 







Table 3. Output morphological type codes 



When the morphological type is uncertain (st > 4.) a 
rough type is used with a question mark (e.g. E? or S?). 
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The histogram of mean morphological type codes is 
given in Fig. 0. These type codes are available for 60130 
galaxies. 



S2 ^ 



Ph o 



5 10 

morphological type code 



Fig. 2. Histogram for morphological type codes. 

The luminosity class introduced by Van den Bergh 
(1960) has been numerically coded between 1 and 9 with 
the extension introduced by H.G. Corwin (SGC, 1985) be- 
tween 9 and 11. The codes are given in Table 5. This code 



code 


original class 


code 


original class 


1 


I 


7 


IV 


2 


III 


8 


IV-V 


3 


II 


9 


V 


4 


II-III 


10 


V-VI 


5 


III 


11 


VI 


6 


III-IV 







Table 4. Luminosity class codes 



is treated as a continuous parameter. We calculated the 
mean luminosity class Ic and the actual uncertainty sic by 
giving the same weight to each reference. 

Luminosity class codes and morphological type codes 
are used to calculate the luminosity index lambda accord- 
ing to de Vaucouleurs et al. (1979). 



lambda = (Ice + t)/10 



(3) 



where Ice is the luminosity class corrected for inclination 
according to the relation: 



Ice = cl — c.logrlh 



(4) 



where logr2h is the axis ratio defined in section 6 and c 
is a parameter depending on the morphological type code 
t (c = 2.0 for t < 5 and c = 2.0 - 0.9 (t - 5), otherwise). 



6. Diameter and axis ratio 

Many papers were devoted to the study of diameters and 
specially to the reduction to the standard system defined 
by the isophote at the brightness of 25B-mag arcsec~ 2 . 
The conclusion of these studies was published by Paturel 
et al. (1991). 

The diameters are expressed to 0.01' in log of 0.1' 
according to the convention of Second Reference Cata- 
log (de Vaucouleurs et al., 1976). They are designated as 
logd2h. For instance a diameter of 10' will be given as 
logd25 = 2.00. Axis ratios are expressed in log of the ratio 
of the major axis to the minor axis. They are designated 
as logr25. 

The main catalogs are reduced to the i?25-standard 
system using a relationship 



logd2b = adogD + b 



logr25 = a' .logR 



(5) 



(6) 



where D is the diameter and R is the ratio of the major 
axis to the minor axis in a given catalog. The constants a, 
b and a' are given in Paturel et al. (1991, tables la and lb) 
for the most common catalogs. Diameters and axis ratios 
extracted from LEDA images or from COSMOS database 
were converted into the standard system using the same 
relationships but with different coefficients (Paturel et al., 
1996; Gamier et al., 1996; Rousseau et al., 1996). 

The completeness curve logN vs. the limiting logd2b 
is shown in Fig. |[ The completeness is satisfied down to 
the limit logdi = 0.9 (i.e. 0.8' in diameter). Diameters 
logd2h are available for 82033 galaxies. The histogram of 
actual uncertainty slogd25 on apparent diameter logd2b is 
given in Fig. ^. More than 13,000 galaxies have a diameter 
with an actual uncertainty smaller than 0.05 (in logd25). 
The distribution of logarithms of axis ratios is shown in 
Fig. H This distribution is close to the one expected if the 
orientation of galaxies is randomly distributed. 

The position angle of the major axis is noted pa. It 
is counted from North towards East, between Odeg and 
180 deg and is almost randomly distributed (Fig. g). A 
small excess of galaxies appears at pa = 90 deg and pa = 
180 deg which seems to be an artifact. 

In RC2 apparent diameters were corrected for galactic 
extinction and inclination effect according to Heidmann, 
Heidmann and de Vaucouleurs (1972abc). Recently, this 
question was revisited after the result by Valentijn (1990, 
1994) that galaxy disks are opaque. Our conclusion (Bot- 
tinclli et al. 1995) leads to the following correction: 



logdc — logd25 — C.logr25 + ag.Ko 



(7) 



where C=0.04, ag is the galactic extinction (see section 
about coordinates) and Ku is given by Fouque and Pa- 
turel (1985) as 0.094 for spiral galaxies and 0.081 - 0.016.i 
for early type galaxies with morphological type code t < 0. 
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1 1.5 2 

limiting diameter logd 1 




0.5 1 

log of axis ratio logr25 



Fig. 3. Completeness curve for logd25. The completeness is 
satisfied down to the limit logdi = 0.9 (i.e. 0.8' in diameter). 



Fig. 5. Histogram of log of axis ratio logr25. The solid curve 
shows the distribution of logr25 for random orientation. 




0.1 0.2 0.3 0.4 

actual uncertainty on logd25 



o ^ - 



r-j - 

p, s - 
o 2 - 



50 100 150 

position angle beta (deg.) 



Fig. 4. Histogram of actual uncertainty slogd25 on apparent 
diameter logd25. 



Fig. 6. Distribution of major axis position angles. 



7. Apparent magnitude and colors 

The reduction of apparent B-magnitude to the RC3- 
system of magnitudes Bt with photoelectric zero-point 
has been studied recently (Paturel et al., 1994). Apparent 
B-magnitudes reduced to the RC3-system will be desig- 
nated as bt. Several effects were taken into account. The 
reduction of a given magnitude to to bt is given by: 



bt = a.m + b+ 

c.(logr25- < logr2b >) 
e.(logd25- < logdlh >) 



d.(t- <t>) + 
-/.((M950- < del950 >) 



(8) 



where a, b, c, d, e, /, < logr25 >,<<>,< logd25 >, 
< del950 > are constant values given in Paturel et al. 
(1994, table 6). The mean bt magnitude is calculated as a 
weighted mean where the weight is derived for each source 
of magnitude as the inverse square of the mean standard 
deviation. The final actual uncertainty sbt is derived from 
the total weight. 

The cumulative completeness curve logN vs. bt is 
shown in Fig. 0. The completeness in apparent magnitude 
is satisfied up to bt = 15.5. Apparent total magnitude bt 
is available for 76760 galaxies. The histogram of actual 
uncertainties sbt on bt is given in Fig. |[ More than 7,000 
galaxies have an actual uncertainty on bt smaller than 0.02 
mag. 
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Note that apparent diameter can be roughly converted 
into a magnitude m assuming that the mean surface 
brightness is constant for all galaxies. The conversion can 
be made using the relationship (Di Nella and Paturel, 
1994): 



20.0 - 5.logd25 



(9) 



The standard deviation on m is about 0.5mag. Using 
this relation, it is possible to obtain an estimate of the ap- 
parent magnitude for 93062 galaxies, 76760 magnitudes of 
which come from bt and 16302 from logd25. This magni- 
tude m will be used for drawing a more general complete- 
ness curve (Fig. |TT|). 



3 
ft 




0.5 



actual uncertainty on bt 



\ x ' " _ Fig. 8. Histogram of actual uncertainty on apparent magni- 

/ | - tudes bt. 

- 

/ (t < 0) we assume ai = 0, in agreement with de Vau- 

\ y / ~_ couleurs et al. (1991). 

- << - Colors are given in the UBV system □. They are: total 

/ ' _ asymptotic colors ubt for (U — B)t ; bvt for (B — V)t and 

x effective colors (i.e. colors within the effective aperture in 

s - which half the total B-flux is emitted) ube for (U — B) e 

>• ; bve for (B — V) e . Total asymptotic colors are corrected 

10 12 14 16 18 for galactic extinction, inclination and redshift effects ac- 
cording to RC3. The corrected colors are bvtc and ubtc for 
(B — V)t and (U — B) T respectively. 



limiting 



Fig. 7. Completeness curve for apparent magnitude bt. The 
completeness is satisfied up to the limit bt = 15.5. 

Apparent bt magnitudes are corrected for galactic ex- 
tinction, inclination and redshift effects according to the 
relation: 

btc = bt-ag-ai- ak.v/WOOO (10) 

where ag is the galactic extinction in B (see section 4), 
expressed in magnitude, ak is taken from de Vaucouleurs 
et al. (1976, RC2 p33, rel.25), v is the heliocentric velocity 
in km.s~ l (section 9) and ai is given by Bottinelli et al. 
(1995) as: 

ai = 2.5log(k + (1 - k).R 2C ^ +0 - 2 ' K ^- 1 ) (11) 

where, k = iBuige/Wotai is taken from Simien and de 
Vaucouleurs (1986), as a function of the morphological 
type code. Kr> is taken from Fouque and Paturel (1985) 
as seen before (section 6), C = 0.04 (Bottinelli et al., 1995) 
and R = 10 a ' r25 . Note that this relation has been demon- 
strated for spiral galaxies only. For early type galaxies 



8. Maximum velocity rotation and central velocity 
dispersion 

We published compilations of Hi-data in 1982 and 1990 
(Bottinelli et al, 1982; Bottinelli et al, 1990) but the data 
are regularly updated from literature. The reduction of 
raw measurements is the same. The 21-cm line widths are 
reduced to two standard levels (20% and 50% of the peak) 
and to zero- velocity resolution using the following formula: 

ws(l; r = 0) = w(l', r) + (af + b)r + c(l - I') (12) 

where ws(l,r = 0) is the standard 21-cm line width at 
level / = 20 or / = 50, while w(l' , r) is a raw measurement 
at a level I' made with a velocity resolution of r km.s . 
The constants a, b and c are (Bottinelli et al., 1990): a = 
0.014,6= -0.83, c= -0.56 

The resulting standard 21-cm line widths ws(l = 
20, r = 0) and ws(l — 50, r = 0) are corrected for sys- 
tematic errors by intercomparison reference by reference 

2 some additional colors V — R, R — I are listed in LEDA 
but they do not represent a significant enough sample to be 
included in the present version of the mean parameters 
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(program INTERCOMP, Bottinelli et al., 1982) leading to 
standard widths w20 and w50 and their actual uncertain- 
ties sw20 and sw50 respectively. 

w20 and u>50 are used to calculate the log of the max- 
imum velocity rotation following the expression. 



logvm =< log(wc) > —log(2sin(incl)) 



(13) 



where incl is the inclination (in degrees) between the po- 
lar axis and the line of sight calculated from the classical 
formula (Hubble 1926): 



sin' 



'(incl) = 



I _ ^Q-2.logr25 
I — \Q~2.logro 



(14) 



0.43 + 0.053i, if -5 < t < 7 (or logro = 
has been obtained from the most flattened 



where logro 
0.38 if t > 7) 
galaxies. 

< log(wc) > is the weighted mean of the logarithm of 
the line widths w20 and w50 corrected for internal velocity 
dispersion. The adopted weight of level 20% is twice the 
weight of level 50% because it is less sensitive to the def- 
inition of the maximum and also because it corresponds 
to larger fraction of the disk. The correction for internal 
velocity dispersion is taken according to Tully and Fouque 
(1985). 



= w 2 + wt 2 (1 - 2e- w2/wr2 ) - 2w.wt(l - e 



2 A 



")(15) 



where w is either w20 or u>50 and wt = 2a z .k(l), assuming 
an isotropic distribution of the non-circular motions o~ z = 
12km. s^ 1 and a nearly Gaussian velocity distribution (i.e. 
fc(20) = 1.96 and fc(50) = 1.13). 

Mean maximum velocity rotation logvm is available 
for 6415 galaxies, from 34,436 individual measurements 
w20 or u;50. 

The actual uncertainty on logvm can be approximated 
by (For the detailed calculation see Bottinelli et al. 1983): 



pi 



o.i 



0.2 



0.3 



0.4 



actual uncertainty slogvm 



Fig. 9. Histogram of the actual uncertainty on maximum ve- 
locity rotation logvm. 



o 



3 



0.05 0.1 
actual error sloes 



0.15 



slogvm = 0.2 



sw 



slogr2h 2 



w 2 ( 10 2/ °S r25 - l) 2 



(16) 



where sw and w are used for (sw20 or sw50) and (w20 or 
w50), respectively. The histogram of slogvm is presented 
in Fig. § 



A preliminary compilation of central velocity disper- 
sions logs was published in 1985 (Davoust et al.,1985) 
and included in our database. This compilation has been 
regularly updated from literature (including compilations 
made by Whitmore et al. 1985, McElroy 1995 and by Prug- 
niel and Simien 1995). Measurements from various refer- 
ences have been homogeneized using the INTERCOMP 
program (Bottinelli et al., 1982). The mean central ve- 
locity dispersion logs is available for 1816 galaxies result- 
ing from 3402 individual measurements. The actual un- 
certainty slogs in log scale is shown in Fig. Uu. 



Fig. 10. Histogram of the actual uncertainty on central veloc- 
ity dispersion logs. 



In Fig. 11 we present the completeness of kinematical 
parameters logvm or logs in comparison with the total 
completeness curve. The completeness is fulfilled up to 
about m = 12.0 mag. 

9. Radial velocities 

Heliocentric radial velocities are obtained from optical or 
radio measurements vopt or vrad, respectively. The orig- 
inal optical compilation was made for the preparation of 
the RC3 catalog (Fouque et al., 1992). Velocities are cor- 
rected for systematic errors from the intercomparison ref- 
erence by reference. The weight is deduced for each refer- 
ence from this comparison. This allows the calculation of 
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limiting magnitude m ] 



Fig. 11. Completeness curve for m. The completeness is sat- 
isfied up to the limit mi = 15.5 (solid line). This limit drop 
to m ~ 14.2 if we impose that the radial velocity is known 
(dashed line), and to m = 12.0 if we impose that either the 
maximum velocity rotation or the central velocity dispersion 
is known (dotted line). 

the actual uncertainty svopt. Radio velocities come essen- 
tially from 21-cm line measurements (plus some additional 
CO measurements). The agreement between different au- 
thors is generally excellent and there is no need of system- 
atic correction. Mean error svrad is given as a function of 
the velocity resolution (Bottinclli et al, 1990). 

From both vopt and vrad we calculate a weighted 
mean heliocentric velocity v, the weights being the inverse 
squares of svopt and svrad respectively. The final weight 
leads to the actual uncertainty sv. When the discrepancy 
between vopt and vrad is larger than lOOOfcm.s -1 , we do 
not calculate the mean heliocentric velocity v but adopt 
instead the velocity having higher weight. Radial velocity 
v is available for 39667 galaxies. 

From this mean heliocentric velocity v we obtain four 
velocities defined with different reference frames. The ve- 
locity corrected to the galactic center vgsr is obtained by 
a correction of the motion of the Sun with respect to the 
local standard of rest (LSR) and a correction of the LSR 
motion with respect to the galactic center. The resulting 
correction is: 

vgsr = v + 232. sin(al2)cos(ab2)+ 

9.cos(al2)cos(ab2) + 7sin(ab2) (17) 

The velocity corrected to the centroid of the Local 
Group vlg has been adopted following Yahil et al (1977): 

vlg = v + 295Asin(al2)cos(ab2)- 

79.lcos(al2)cos(ab2) - 37.6sin(ab2) (18) 



This correction replaces the classical IAU correction 
300sin(al2)cos(al2). 

The velocity corrected for infall of the Local Group 
towards Virgo is noted vvir. It is calculated as: 

vvir = vlg + I70.cos(9) (19) 

where 170km. s^ 1 is the infall velocity of the Local 
Group according to Sandage and Tammann (1990) and 
where 9 is the angular distance between the observed di- 
rection sgl, sgb in supergalactic coordinates and the di- 
rection of the center of the Virgo cluster (sglo = 104 deg, 
sgbo — — 2 deg). 

cos{9) — sin(sgbo)sin(sgb)+ 

cos(sgbo)cos(sgb)cos(sglo — sgl) (20) 

Finally, the radial velocity is also expressed in the ref- 
erence frame of the Cosmic Background Radiation. This 
velocity is noted v3k. It is calculated from the heliocen- 
tric velocity v using the total solar motion of 360km. s^ 1 
towards the direction defined by the 1950- equatorial co- 
ordinates al3k = 11.25/i de3k = —5.6 deg (Lubin and Vil- 
lela, 1986). In 1997 this calculation should be replaced by 
the new determination from COBE (Bennett et al. 1996). 
However, according to the rule defined at the end of the 
present paper (see the section "Acknowledgements"), the 
old definition will be used until the end of 1996: 

v3k = v + 360.cos(93k) (21) 

with 93k given by: 

cos(93k) = sin(de3k)sin(del950) + 

cos(de3k)cos(del950)cos(al3k - aZ1950) (22) 

10. Hl-line and IR fluxes 

HI line flux (flux corresponding to the area under the 21- 
cm line profile) and IRAS fluxes are treated separately 
from the classical magnitudes (UBV) because they are 
obtained and corrected in a completely different way. 

The HI line flux is generally expressed in Jy.km.s~ 1 
converted in magnitude m21 according to the formula 
adopted in RC3: 

m21 = -2.5log(f) + 17 .40 (23) 

where / is the area of the 21-cm line profile expressed 
in Jansky.km.s~ 1 . This formula is equivalent to the one 
used in RC3: m21 = —2.5log(fwm) + 16.6, where fwm 
is the flux in 10~ 22 W.m~ 2 . The standard error sm21 on 
m21 is given as a function of the radiotelescope according 
to Bottinelli et al. (1990). 
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The HI line magnitude m21 has been corrected for self- 
absorption effect following Heidmann et al. (1972): 

K 1 1 cos{incl) 



m2\c — m21 — 2.5log- 



(24) 



(1 — exp(/«/ cos{incl)) 

The adopted free parameter is n = 0.031. If inclination 
is higher than 89deg the maximum correction is limited 
to — 0.82mag. 

Both magnitudes m21c and btc are used to calculate 
a HI color index hi initially defined by de Vaucouleurs et 
al. (1976): 



hi ~ to21c — btc 



(25) 



This parameter is interesting as it is directly connected 
to the hydrogen contents per unit of B-flux (see RC3, p51 
Rel. 78). The relation between hi and morphological type 
code t is presented in Fig. |l2|. It shows a clear correlation 
which validates the use of morphological type code as an 
observable parameter. 



aperture enclosing one-half the total flux) . This mean sur- 
face brightness is expressed in B-mag .arcsec.~ 2 . 

Two equivalent measurements of brief were derived: i) 
from the apparent diameter D n enclosing a mean surface 
brightness of20.75B-mag.arcsec.~ 2 (Dressier et al. 1987), 
ii) from the mean surface brightness inside the effective 
isophote (elliptical isophote enclosing one-half the total 
flux) measured by Lauberts and Valcntijn (1989; LV). 

brief = -(1.87 ± 0.07)m'(-D„) + (60.13 ± 1.36) 

a = 0.49 

p = 0.74 ±0.02 

n = 392 

(27) 

where m'(D n ) = bt+5log(D n )+4:.38, a is the standard 
deviation, p is the correlation coefficient and n the number 
of galaxies used for the comparison. 

Similarly we have: 



- -O 



T3 



-o 



- (X 



10 



morphological type code t 



Fig. 12. Mean HI color index hi vesus morphological type code 
t. 



IRAS fluxes at QQp,m and 100/xm are converted in the 
so-called far-infrared flux according to the relation: 



mfir = -2.5; O <?(2.58./60+ /100) + 14.75 



(26) 



where /60 and /100 are IRAS fluxes at 60/im and 
100/im expressed in Jansky. This relation is equivalent to 
the relation given in RC3. The term 14.75 comes from 
the arbitrary zero-point of 20 in RC3 (p.43, Rel. 49). The 
factor 2.58 comes from IRAS Point Source Catalog (1988). 



brief = (1.14 ± 0.02)m'{LV) + (3.05 ± 0.46)+ 

(1.2±0.1)Zogr25 2 

a = 0.49 

p = 0.84 ±0.01 

n = 847 

(28) 

where m'(LV) is the average blue central surface 
brightness within half total B light (noted p e (LV) in LV) 
, from Lauberts and Valentijn (1989). The correction for 
inclination is in good agreement with the predicted one 
1.35 - 1.26 (see RC3 p50, Rel.71). The final value brief 
is calculated as the weighted mean of each determination. 

Another estimate of the mean surface brightness is 
6ri25, the mean surface brightness inside the isophote 
2hB-mag .arcsecr 2 . 

This brightness must be corrected for inclination ef- 
fect (Bottinclli et al., 1995). We then obtain the corrected 
mean brightness bri2b: 

bri25 = to'25 + 2.5log(k.R- 2C + (1 - k)R^ OAC ^ K ^- 1 ){29) 
where 



i'25 = bt + 5 logd25 + 3.63 



(30) 



Notations are those used for the calculation of ai (sec- 
tion 7). 



11. Mean surface brightness 

The parameter brief (with its actual uncertainty sbrief) 
is the mean effective surface brightness, i.e. the mean sur- 
face brightness inside the effective aperture (the circular 



12. Distance modulus and absolute magnitude 

The kinematical distance modulus mucin can be derived 
from the heliocentric velocity properly corrected for the 
Local Group infall onto the Virgo cluster vvir assuming 
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a given Hubble constant H = 75km. s^ 1 .Mpc^ 1 . It must 
be noted that this distance does not include a correction 
for the infall of individual galaxies onto Virgo. Such a dis- 
tance could have been calculated for instance using the 
model by Peebles (1976) as described in Bottinelli ct al. 
(1986). However, this model does not allow the calculation 
in the direction of the Virgo center because of the third 
degree equation (Eq. 2 in Bottinelli et al. 1986). Further, 
this model requires the choice of a Virgo distance, of a 
Virgo mean radial velocity and of a velocity infall for the 
Local Group, while the calculation of mucin requires only 
the choice of the velocity infall for the Local Group (we 
adopted Vinfaii = 170km. s^ 1 ; see section 9). In the direc- 
tion of Virgo cluster center mucin can be overestimated 
or underestimated depending on the background or fore- 
ground position of the considered galaxy with respect to 
Virgo center. 

mucin — 5.log(vvir/75) + 25 (31) 

mucin is calculated only where vvir > 500km. s^ 1 , 
mucin is available for 39243 galaxies. It is used to de- 
rive an estimate of the absolute magnitude amabs in Blue 
band: 

mabs = btc — mucin (32) 
13. Group membership 

The extragalactic database was used by Garcia (Garcia et 
al. 1993; Garcia 1993) for a general study of group mem- 
bership of all galaxies with a radial velocity less or equal 
to 5500 km.s' 1 and an apparent magnitude brighter than 
bt = 14. Two automatic algorithms were used simultane- 
ously (percolation and hierarchy clustering methods) for 
producing very robust groups. From the whole sample 485 
groups were build. They arc identified by the acronym 
LGG (for Lyons Galaxy Group). For a galaxy, the LGG 
number Igg gives the group to which the galaxy belongs. 
The Igg number is available for 2702 galaxies. 
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in the LEDA news. The data of previous years will be 
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Appendix A: Calculation of the actual uncertainty. 

Let us assume that for a given galaxy we have n measure- 
ments Xi (i = 1, n) of a given parameter obtained from different 
references, each reference having a weight u>; = 1/of , where Ui 
is the standard error of the i-th individual measurement. The 
actual uncertainty is calculated as: 

Ca.u. — &w + °"n (1) 

The first term a w denotes the inverse of the total weight 
The total weight, S w , is simply the sum of individual weights, 
(i.e. S w — l/cr w — J^l/of). This first term accounts for the 
accuracy of the reference of each individual measurement be- 
cause the standard error of a given measurement is assigned 
globaly from e.g., the reference or the resolution etc... It is ob- 
vious that some individual measurements coming from a good 
reference can be affected by a local problem (e.g., multiplic- 
ity of the galaxy, star superimposed on the galaxy, bad seeing, 
misidentification etc.). This fact will be taken into account by 
the second term. 

The second term in the definition of the actual error is a 
measure of the consistency of the different measurements build- 
ing the mean measurement. It is calculated as the weighted 
standard deviation: 

On = S2/S w — (Sl/S w ) 2 (2) 

where : 5 2 = £V Wix\, Si = £\ Sw = J2i Wi - 

The main advantage of the actual error is that it clearly 
shows any internal uncertainty and any external discrepancy. 
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Appendix B: LEDA's Astrophysical parameters 



Parameter FORTRAN FORTRAN columns definition 
name name format 



pg c 


pgcleda 


all 




1- 


- 11 


PGC or LEDA name (PGC = LEDA) 


ident 


ident 


al6 




12- 


- 27 


1st name (NGC , IC , UGC , ES0 . . . ) 


ipad 


ipadc 


al 




28- 


- 28 


'*' for coordinates better than 10" 


all950 


all950 


flO. 


.5 


29- 


- 38 


R.A. (B1950) (decimal hours) 


del950 


del950 


f 10. 


.5 


39- 


- 48 


DEC. (B1950) (decimal degrees) 


al2000 


al2000 


f 10. 


.5 


49- 


- 58 


R.A. (J2000) (decimal hours) 


de2000 


de2000 


f 10. 


.5 


59- 


- 68 


DEC. (J2000) (decimal degrees) 


12 


al2 


flO. 


.3 


69- 


- 78 


galactic longitude (degrees) 


b2 


ab2 


flO. 


.3 


79- 


- 88 


galactic latitude (degrees) 


sgl 


sgl 


f 10. 


3 


89- 


- 98 


supergalactic longitude (degrees) 


sgb 


sgb 


flO. 


.3 


99- 


-108 


supergalactic latitude (degrees) 


"typ 


type 


a5 , lx 


109- 


-114 


morph. type (e.g. 'E' , 'Sab' , 'SBa' , 'SO' ) 


morph 


morphc 


a4 




115- 


-118 


'B' for Barred gal. (see note below) 














'R' for Ring gal. 














'M' for multiple gal. 














'C for compact, 'D' for diffuse 


t 


t 


f 10. 


3 


119- 


-128 


morph. type code (-5 to 10) 


St 


St 


f 10. 


.3 


129- 


-138 


actual uncertainty on t 


lc 


ale 


flO. 


.3 


139- 


-148 


luminosity class (1 to 11) 


sic 


sic 


f 10. 


.3 


149- 


-158 


actual uncertainty on lc 


logd25 


alogd25 


f 10. 


.3 


159- 


-168 


loglO of isophotal diameter (d25 in 0.1') 


slogd25 


slogd25 


flO. 


.3 


169- 


-178 


actual uncertainty on logd25 


logr25 


alogr25 


f 10. 


.3 


179- 


-188 


loglO of the axis ratio (major/minor axis) 


slogr25 


slogr25 


f 10. 


3 


189- 


-198 


actual uncertainty on logr25 


pa 


pa 


f 10. 


.3 


199- 


-208 


position angle (N->E) in degrees 


brief 


brief 


f 10. 


.3 


209- 


-218 


effective surface brightness (mag. arcsec-2) 


sbrief 


sbrief 


flO. 


3 


219- 


-228 


actual uncertainty on brief 


bt 


bt 


f 10. 


.3 


229- 


-238 


total B-magnitude 


sbt 


sbt 


f 10. 


.3 


239- 


-248 


actual uncertainty on bt 


ubt 


ubt 


f 10. 


.3 


249- 


-258 


(U-B)T 


bvt 


bvt 


flO. 


3 


259- 


-268 


(B-V)T 


ube 


ube 


flO. 


.3 


269- 


-278 


(U-B)e 


bve 


bve 


f 10. 


.3 


279- 


-288 


(B-V)e 


w20 


w20 


f 10. 


3 


289- 


-298 


21-cm line width at 20'/, of peak (in km/s) 


sw20 


sw20 


f 10. 


.3 


299- 


-308 


actual uncertainty on w20 


w50 


w50 


flO. 


.3 


309- 


-318 


21-cm line width at 50'/, of peak (in km/s) 


sw50 


sw50 


flO. 


.3 


319- 


-328 


actual uncertainty on w50 


logs 


alogs 


f 10. 


.3 


329- 


-338 


log of the central velocity disp. (s in km/s 


slogs 


slogs 


flO. 


.3 


339- 


-348 


actual uncertainty on logs 


m21 


am21 


flO. 


.3 


349- 


-358 


Hi-magnitude 


sm21 


sm21 


flO. 


.3 


359- 


-368 


actual uncertainty on m21 


mf ir 


amf ir 


flO. 


.3 


369- 


-378 


far-infrared magnitude 


vrad 


vrad 


f 10. 


.3 


379- 


-388 


radio heliocentric radial velocity in km/s 


svrad 


svrad 


flO. 


3 


389- 


-398 


actual uncertainty on vrad 


vopt 


vopt 


flO. 


.3 


399- 


-408 


optical heliocentric radial velocity in km/: 


svopt 


svopt 


flO. 


.3 


409- 


-418 


actual uncertainty on vopt 


V 


V 


flO. 


.3 


419- 


-428 


actual heliocentric radial velocity in km/s 


sv 


sv 


flO. 


3 


429- 


-438 


actual uncertainty on v 
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Parameter FORTRAN 


FORTRAN columns 


definition 


name name 


format 





!gg 


algg 


flO. 


.3 


439- 


-448 


Lyon's galaxy group number 




ag 


ag 


flO. 


.3 


449- 


-458 


galactic extinction in B-mag 




ai 


ai 


f 10. 


.3 


459- 


-468 


internal absorption (in B-mag) 




incl 


aincl 


f 10. 


.3 


469- 


-478 


inclination 




a21 


a21 


f 10. 


3 


479- 


-488 


HI self -absorption 




lambda 


alambda 


flO. 


.3 


489- 


-498 


luminosity-index 




logdc 


alogdc 


flO. 


3 


499- 


-508 


log of the corrected diameter (dc in 0.1') 


btc 


btc 


f 10. 


.3 


509- 


-518 


corrected B-magnitude 




ubtc 


ubtc 


f 10. 


.3 


519- 


-528 


(U-B)o 




bvtc 


bvtc 


f 10. 


3 


529- 


-538 


(B-V)o 




bri25 


bri25 


flO. 


.3 


539- 


-548 


mean surf . brightness within 25 


m/» 


logvm 


alogvm 


flO. 


.3 


549- 


-558 


log of max.circ. rot. vel . 




slogvm 


slogvm 


f 10. 


.3 


559- 


-568 


actual uncertainty on logvm 




m21c 


am21c 


flO. 


.3 


569- 


-578 


corrected Hi-magnitude 




hie 


hie 


flO. 


3 


579- 


-588 


HI color index 




vlg 


vlg 


flO. 


.3 


589- 


-598 


radial vel . relative to the LG 




vgsr 


vgsr 


flO. 


.3 


599- 


-608 


radial vel. relative to the GSR 




vvir 


vvir 


f 10. 


.3 


609- 


-618 


radial vel. corrected for Virgo 


infall 


v3k 


v3k 


flO. 


.3 


619- 


-628 


radial vel . relative to the CBR 




mucin 


amucin 


flO. 


.3 


629- 


-638 


kinematical distance modulus (H= 


=75 km/s/Mpc) 


mabs 


amabs 


flO. 


.3 


639- 


-648 


absolute B magnitude from mucin 


and mupar 


identi 


identi 


20al6 


649- 


-968 


alternate names 





note: The parameter 'morph' can be read as 4 parameters (4al format) 
for Bar, Ring, Multiple and Compactness, respectively 



This article was processed by the author using Springer- Verlag M?gX A&A style file L-AA version 



